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Abstract 

The low-temperature properties of systems characterized by a sponta¬ 
neously broken internal rotation symmetry, 0(A) —)■ O(A-l), are governed 
by Goldstone bosons and can be derived systematically within effective La- 
grangian held theory. In the present study we consider systems living in 
two spatial dimensions, and evaluate their partition function at low tem¬ 
peratures up to three-loop order. Although our results are valid for any 
such system, here we use magnetic terminology, i.e., we refer to quantum 
spin systems. We discuss the sign of the Goldstone boson interaction in 
the pressure, staggered magnetization, and susceptibility as a function of 
an external staggered held for general N. As it turns out, the d=2+l 
quantum XY model (A=2) and the d=2+l Heisenberg antiferromagnet 
(A=3), are rather special, as they represent the only cases where the 
spin-wave interaction in the pressure is repulsive in the whole parameter 
regime where the ehective expansion applies. Remarkably, the d=2+l XY 
model is the only system where the interaction contribution in the stag¬ 
gered magnetization (susceptibility) tends to positive (negative) values at 
low temperatures and weak external held. 
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1 Introduction 


The present study is devoted to (pseudo-)Lorentz-invariant systems that are defined 
in two spatial dimensions and - at T=0 - are characterized by a spontaneously broken 
rotation symmetry 0{N) —)■ O(A^-l). We focus on the low-energy dynamics which is 
dominated by the corresponding iV-1 Goldstone bosons. The most prominent physical 
realizations are the d=2+l quantum XY model (iV=2) and the d=2+l Heisenberg 
antiferromagnet {N=3). In either case the relevant low-energy degrees of freedom are 
the spin waves or magnons. 

Although there exists a vast number of publications - in particular for iV={2,3} ~ 
the impact of the Goldstone boson interaction onto the low-temperature behavior of 
these systems has been largely neglected. While the specific cases Y={2,3} have been 
addressed within effective Lagrangian field theory in Refs. [H-lsI, a systematic study of 
the nature of the interaction in systems with general N, seems to be lacking. It is the 
aim of the present investigation to rigorously analyze the effect of the Goldstone boson 
interaction in the pressure, in the staggered magnetization, and in the susceptibility, 
using the systematic and model-independent effective Lagrangian technique. 


A method widely used to analyze the low-temperature behavior of spin systems 
with a spontaneously broken internal rotation symmet^ _^spin-wave theory, or its 

With respect to the 


modified versions adapted to two spatial dimensions ^ _ 

d=2+l quantum XY model, the interested reader may also consult Refs. [li-H 
The finite-temperature properties of the d=2+l Heisenberg antiferromagnet and the 
d=2+l quantum XY model, have furthermore been analyzed with alternative methods 
that include Schwinger-boson mean-field theory. Green’s functions, renormalization 
group methods, and Monte-Garlo simulations. Goncerning the d=2+l quantum XY 


model we refer to Refs. |23N33|. while for the d=2+l Heisenberg antiferromagnet we 


point to Refs. j34l-l47l|. For the general case, i.e., for arbitrary iV, see Refs. j48l-l56 


In the work reported here, we do not make use of spin-wave theory or any other mi¬ 
croscopic method. Rather, our approach relies on the effective Lagrangian technique 
which systematically captures the low-energy physics of the system that is governed 
by the Goldstone bosons. Low energy - or low temperature - means that we are 
interested in the properties of the system below the intrinsic scale defined by the 
underlying microscopic model. In connection with the Heisenberg antiferromagnet or 
the quantum XY model, this scale is given by the exchange coupling J. Note that 
the spontaneously broken rotation symmetry 0(iV) is only approximate, since we in¬ 
corporate a weak external field into the low-energy description. As we discuss below, 
this field cannot be completely switched off in two spatial dimensions in the effective 
approach we pursue in the present study. 

The d=2+l quantum XY model {N=2) and the d=2+l Heisenberg antiferromag¬ 
net (Y=3) turn out to be rather peculiar cases: they represent the only systems where 
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the spin-wave-interaction contribution in the pressure is repulsive at low temperatures 
in the entire parameter regime where our effective analysis applies. Furthermore, con¬ 
sidering the impact of the spin-wave interaction in the staggered magnetization and 
susceptibility, the low-temperature behavior of the d=2+l quantum XY model is 
quite different from any other d=2+l (pseudo-)Lorentz-invariant system with a spon¬ 
taneously broken internal rotation symmetry 0{N): it is the only system where the 
interaction contribution in the staggered magnetization (susceptibility) tends to pos¬ 
itive (negative) values at low temperatures and weak external field. Interestingly, for 
systems with N > 4 the staggered magnetization (susceptibility) also becomes posi¬ 
tive (negative) at very low temperatures and not too weak external helds. This means 
that here the Goldstone boson interaction at hnite temperature tends to reinforce the 
antialignment of the spins and enhances the staggered spin pattern. 

We emphasize that the physical realizations for Y={2,3} are not restricted to 
quantum magnetism. Rather, any two-dimensional system with a spontaneously bro¬ 
ken rotation symmetry, 0(3) —)■ 0(2), or 0(2) —)■ 1, exhibits the same characteristic 
features determined by the universal physics of its Goldstone bosons. Most remark¬ 
ably, the question of whether the Goldstone boson interaction in the pressure and the 
other observables we consider, is repulsive or attractive at low temperatures, applies to 
any d=2+l (pseudo-)Lorentz-invariant system that exhibits a spontaneously broken 
rotation symmetry. In this sense, our results on the nature of the Goldstone boson 
interaction are universal. To the best of our knowledge, our three-loop results are 
far beyond any analogous calculation based on more conventional condensed matter 
methods like spin-wave theory. 


Regarding the d=2+l quantum XY model and the Heisenberg antiferromagnet, 
high-precision measurements based on very efficient loop-cluster algorithms have been 
performed recently in Refs. [^. These numerical results were confronted with 


analytic calculations obtained within effective Lagrangian theory. Not only were some 
combinations of low-energy constants determined at the permille level, but also was 
the correctness of the systematic effective Lagrangian method demonstrated. 


In order to convince the interested reader that the effective Lagrangian technique 
indeed represents a very powerful machinery in the condensed matter domain, we 
point out the following additional references which are all based on that method. 
Ferromagnets and antiferromagnets in one, two, and three spatial dimensions have 
been analyzed in Refs. (m-^. (Quasi-)two-dimensional antiferromagnets that turn 
into high-temperature superconductors upon doping with either holes or electrons, 
have been systematically studied in Refs. (79-90|. 


The rest of the paper is organized as follows. In Sec. [2] we give a concise account 
of the effective Lagrangian method. The three-loop analysis of the partition function 
is reviewed in Sec. [3l The low-temperature series for the pressure, the staggered 
magnetization, and the susceptibility for d=2+l (pseudo-)Lorentz-invariant systems 
with a spontaneously broken symmetry, 0(iV) —)■ O(iV-l), are discussed in Sec. 01 Our 
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focus is the manifestation of the Goldstone boson interaction in the low-temperature 
behavior of these systems in a weak external held. Finally, in Sec. Elwe present our 
conclusions. 


2 Effective Field Theory 


The basic principles underlying the effective Lagrangian method, as well as the per¬ 
turbative evaluation of the partition function, have been outlined on various occasions 
and will not be repeated here in detail. Rather, in this section, we discuss some es¬ 
sential aspects of the method, such that a non-expert can understand the calculation 


presented here. The interested reader may hnd details in Appendix A of Ref. 


or in Section 2 of Ref. [TOj]. Moreover, for pedagogic introductions to the effective 


Lagrangian technique we refer to Refs. j9ll-l96 


The basic degrees of freedom in the effective Lagrangian are the N-1 Goldstone 
bosons, which result from the spontaneously broken rotation symmetry 0{N) 
0{N-1). These Goldstone boson fields are collected in a V-dimensional unit vector 
W{x) 

U\x) U\x) = 1, i = 1, 2,..., V. (2.1) 

While the hrst N-1 components describe the Goldstone bosons, the last component, 

U^(^x) = y/l-G“(a;)f/“(a:), a = 1, 2,..., V - 1, (2.2) 


guarantees that the vector U = 17* (x) is normalized to unity. 

The low-energy dynamics of the system, dominated by the Goldstone bosons, is 
described by the following Euclidean effective Lagrangian: 

= C-lff + , 

(2.3) 

The different pieces in Ceff are organized according to the number of derivatives they 
contain. The leading contribution L^e// contains two space-time derivatives (order 
p^), while the next-to-leading contribution is of order p^. The external field if* 
counts as a quantity of order p^. 

A physical realization for N=2 and N=3 is provided by the 7=2-1-1 quantum XY 
model and the 7=2-|-l Heisenberg antiferromagnet, respectively. The 0(2) spin sym¬ 
metry of the quantum XY Hamiltonian, and the 0(3) spin symmetry of the Heisenberg 
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Hamiltonian, are spontaneously broken according to 0(2) —)■ 1 and 0(3) —)■ 0(2), re¬ 
spectively. The Goldstone bosons in these cases are identihed with the spin waves or 
magnons, and the external held Hs = ( 0 ,..., 0 , Hg) is the staggered held that couples 
to the staggered magnetization vector U. Note that there exists a mapping between 
the antiferromagnetic and the ferromagnetic quantum XY model on a bipartite lat¬ 
tice, such that for N=2 the external held may also be interpreted as the magnetic 
held that couples to the magnetization vector U of the quantum XY ferromagnet. 

It should be stressed that the ehective description is universal, in the sense that 
any two (pseudo-)Lorentz-invariant systems that display the same internal and space- 
time symmetries, will be described by the same ehective Lagrangian. But while the 
structure of the ehective Lagrangian is the same, the diherence between two distinct 
physical realizations manifests itself in the actual numerical values of the low-energy 
coupling constants F, ci, 62 , fci, ^ 2 , ^ 3 - Symmetry leaves these ehective constants 
undetermined, such that they have to be extracted from a numerical simulation, by 
comparison with the microscopic calculation, or by experiment. In the context of 
quantum spin systems, is referred to as helicity modulus or spin stihness, while 
Ss is the staggered magnetization at zero temperature. 

An important point to emphasize is that the ehective Lagrangian, Eq. fl2.3jl . is 
written in a Lorentz-invariant form. One may thus be very skeptical about its appli¬ 
cability to the d=2+l quantum XY model or the Heisenberg antiferromagnet, whose 
underlying two-dimensional lattices are not even invariant under continuous space ro¬ 
tations. First of all, the leading piece does not depend on the specihc geometry 
- the lattice structure only starts manifesting itself at subleading orders in the ehec¬ 
tive Lagrangian j^. Hence may be written in a (pseudo)-Lorentz-invariant form, 
where the spin-wave velocity takes over the role of the velocity of light. Note that 
(pseudo-)Lorentz-invariance is an accidental symmetry of Fg^^:, and is not present in 
the microscopic XY or Heisenberg Hamiltonian. As far as the next-to-leading piece 
Fgjj is concerned, indeed, lattice anisotropies start to show up at this order: these 
contributions are still invariant under the discrete spatial symmetries of the lattice, 
but break continuous space rotation invariance and thus (pseudo-)Lorentz invariance. 
As will become clear in the next section, these extra contributions do not ahect our 
conclusions concerning the impact of the spin-wave interaction in thermodynamic 
quantities. In other words, it is perfectly safe for our purposes to also write Fg^j in a 
(pseudo-)Lorentz-invariant form. 

Since we are dealing with relativistic kinematics, the external held Hg induces a 
mass term - or energy gap - in the dispersion relation of the Goldstone bosons, 

a;= + = (2.4) 

where v is the spin-wave velocity. The low-energy degrees of freedom are then re¬ 
ferred to as pseudo-Goldstone bosons. As long as the external held Hg, and therefore 
M^, are small compared to the intrinsic scale occurring in the underlying micro- 
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scopic model, these quantities may be treated as perturbations within the effective 
Lagrangian framework. In the case of the quantum XY model or the Heisenberg 
antiferromagnet, this scale is given by the exchange integral J that appears in the 
corresponding microscopic Hamiltonians. 


3 Three-Loop Analysis of the Partition Function 


The evaluation of the partition function for (pseudo-)Lorentz-invariant systems, de¬ 
fined in two spatial dimensions and exhibiting a spontaneously broken symmetry 0(iV) 
—)■ O(Y-l), has been presented in Ref. within effective Lagrangian field theory up 
to the three-loop level. In order not to be repetitive, here we restrict ourselves to a 
minimum of information - just to understand the subsequent discussion. 


Up to three-loop order, one has to evaluate the eight Feynman diagrams depicted 
in FigJU We should stress that we are dealing with a systematic loop expansion 
of the partition function, where each Goldstone boson loop is suppressed by one 
power of momentum. The two-loop graph 4b, compared to the one-loop graph 3, is 
suppressed by one additional power of p. The three-loop graphs 5a-c, on the other 
hand, are suppressed by two additional powers of momentum compared to graph 3, 
such that they are of next-to-next-to leading order p^. While the evaluation of most of 
these graphs turns out to be straightforward within effective field theory, the piece of 
resistance is the three-loop graph 5c. The renormalization and numerical evaluation 
of this graph has been presented in detail in Refs. [3, • 

The final result for the free energy density {N > 2) up to three-loop order is 


2 = 2„-i(iv-i)V<2)r’+ 

- j^(N-l)(N-3f^hdafh,(a}T^ + -^q(a)T’’ + 0{T<‘), (3.1) 

where the ratios 

T 


define dimensionless variables, and zq is the free 


Mn 2 ^ 

energy density at zero temperature. 
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Figure 1: Feynman diagrams concerning the low-temperature expansion of the par¬ 
tition function for systems with a spontaneously broken rotation symmetry 0{N) —)■ 
O(A^-l), up to three-loop order and for two spatial dimensions. The number 4 at¬ 
tached to the vertices in graphs 4a and 5d, indicates that they originate from 
Vertices involving the leading piece are denoted by a hlled circle. In d=2+l, 
each Goldstone boson loop is suppressed by one momentum (or temperature) power. 


The quantity is the renormalized mass of the N-1 pseudo-Goldstone bosons, 

EsHs {N - 3) 


= 


F2 


Svr 

^i) + 


F5 


hi 


F2 + 647r2F2 / F6 


.1 


S2V2 




(3.3) 


or, equivalently, the energy gap in the presence of the external held Hg. Note that 
the hrst term on the RHS is the leading mass contribution (order oc Hg) dehned in 
Eq. fl2.4p . while the remaining two contributions in Eq. fl3.3p are subleading corrections 
of order and respectively. The quantities ki and /c 2 are the next-to-leading-order 
effective constants dehned in Eq. fl2.3p . and the explicit expressions for the coefficients 
hi and 62 are listed in Appendix B of Ref. j^. 


The dimensionless kinematical functions hg, hi and h 2 appearing in Eq. fl3.ip . are 
proportional to Bose functions gj.{M,T) in d space-time dimensions, 

go = T^ho{a), gi = T hi{a), g^ = ^ h2{a), (3.4) 

which are dehned as 

/»00 J 

gr{M,T) = 2 exp(-pM2) ^ exp(-uV4pT2), (3.5) 

(47rp) 
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Figure 2: The kinematical functions hk{cr)\ k=Q (continuous curve), k=l (dashed 
curve), and k=2 (dotted curve), in terms of the parameter a = M^/2 ttT. 


and which we have depicted in Fig. [2J Note that in the limit a —)■ 0, the functions hi 
and ^2 diverge, whereas ho tends to a hnite value. 

Finally, the function q appearing in the last term of Eq. 03.1 p is given by 

= 5^(iV-l)(iV-3)M*J,-i(JV-l)(iV-2) J,. 

= 5k(iV-l)(iV-3)<;i(<T)r=-i(]V-l)(]V-2)9jMr=. (3.6) 

A plot for the two individual functions qi and q 2 , in terms of the parameter a = 
Mt^/2ttT^ is provided in Fig. |3l The function q involves the renormalized integrals Ji 
and J 2 , 


Ji = Ji - Cl - 02 5-1 (M,T), 

J2 = <^2 — C3 — C4 5-i(M, T), ( 3 . 7 ) 


stemming from the three-loop graph 5c. These hnite integrals are related to the 
unrenormalized integrals Ji and J 2 over the torus T = TZ'^~ ^ X via 


Ji 

J2 




5 




'T 


df,G{x) d^G{x] 
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(3.8) 




qiio-) q\{o-) 




Figure 3: The three-loop functions q'i(cr) (top) and q' 2 (cr) (bottom), in terms of the 
parameter a = M^/27rT. 
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The coefficients ci,..., C 4 represent subtraction constants which remove the infinities 
in the integrals Ji and J 2 - for the explicit expressions the reader may consult Refs. [Tl, 

a. 


The low-temperature series for the free energy density of (pseudo-)Lorentz-invariant 
systems, defined in two spatial dimensions and displaying a spontaneously broken sym¬ 
metry O(A^) —)■ O(A^-l), thus exhibits the following general structure. The leading 
contribution is of order - this is the free Bose gas term which originates from the 
one-loop graph 3. The interaction sets in at order and corresponds to the two-loop 
graph 4b. We then have four interaction contributions of order T® that come from 
graphs 5a-c. 

Note that the one-loop graph 5d that describes noninteracting Goldstone bosons, 
is the only loop graph that involves an insertion from where lattice anisotropies 
start to show up - it is important to stress that the low-temperature structure of 
the interaction terms is therefore not affected by these anisotropies. There exists, 
however, a subtle effect as the one-loop graph 5d manifests itself in the renormalized 
Goldstone boson mass given in Eq. fl3.3jl : the combination ^2 - ki of next-to-leading 
effective constants in the curly bracket originates from graph 5d, indicating that lattice 
anisotropies start showing up at this order. But note that we are dealing with a next- 
to-next-to-leading order effect in a temperature-mdependent quantity which does not 
affect our conclusions regarding the nature of the Goldstone boson interaction. 


4 Nature of Goldstone Boson Interaction 


Let us explain what we mean by nature of the Goldstone boson interaction. First 
of all, we refer to a diagrammatic definition of the interaction. Tree graphs or one- 
loop graphs related to the partition function describe free Goldstone bosons. But any 
partition function diagram that contains two or more loops, represents an interaction 
diagram: here, four or more Goldstone boson lines meet at one vertex. In the present 
context of d=2+l (pseudo-)Lorentz-invariant systems with a spontaneously broken 
rotation symmetry, the interaction manifests itself at two-loop order (diagram 4b) 
and at three loop-order (diagrams 5a-c). 

Then we adopt a thermodynamically motivated interpretation of the interaction, 
which refers to the sign of the interaction contribution in the low-temperature series of 
the various quantities we consider - this is what we mean by nature of Goldstone boson 
interaction. Regarding the pressure, a positive (negative) interaction contribution is 
interpreted as repulsive (attractive), since it is the picture of the non-ideal Goldstone 
boson gas we have in mind. Goncerning the low-temperature series for the staggered 
magnetization (order parameter) and the staggered susceptibility, the sign of the 
interaction contribution is related to the antialignment of the spins, as we will discuss. 
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4.1 Effective Expansion: Domain of Validity 


Before submerging into the analysis of the interaction, we should point out in which 
parameter regime - dehned by T, Hg (or M^-), F, J - the effective series are valid. 
To explore this domain, we introduce the two dimensionless ratios 


T Hg 

2ttF^ ’ 27rF2 ' 


(4.1) 


Note that, at leading order, the external held Hg is proportional to the square of the 
Goldstone boson mass according to Eq. fl3.3p . 


For the effective expansion to be valid, both the temperature and the external 
held must be small with respect to the underlying scale of the system. In chiral per¬ 
turbation theory, i.e., the ehective held theory of quantum chromodynamics (QCD), 
this underlying scale - the chiral symmetry breaking scale Aqcd ^ is of the order of 
Aqcd = where F ~ 93MeV is the pion decay constant. In the present study, 

we adopt this dehnition of the scale, bearing in mind that in connection with spin 
systems with a spontaneously broken rotation symmetry (and in two spatial dimen¬ 
sions), the analog of the pion decay constant is the helicity modulus or spin stihness 
This quantity has been measured in high-precision Monte Carlo simulations, both 
for the d=2+l quantum XY model and the d=2-|-l Heisenberg antiferromagnet on 
the square lattice, with the result 26, 0, 58|: 


= 0.26974(5) J, = 0.1808(4) J. (4.2) 


Here the underlying scale - the exchange constant J of the microscopic theory - can 
thus be expressed in terms of the effective coupling constant F as J ~ 2ttF‘^ [^. 
Accordingly, the two dimensionless quantities dehned in Eq. fl4.ip simply are 


T 

7 ’ 



(4.3) 


In order to determine what ’’low temperature” means, let us again consider the 
situation in quantum chromodynamics. Unlike in the present case where <7=2+1, in 
QCD we are dealing with three spatial dimensions. QCD is characterized by a chiral 
phase transition taking place at Tqcd ~ lOOMeU where the Coldstone boson picture 
dehnitely ceases to be valid. Compared to the scale of QCD, Aqcd ~ IGeU, we 
obtain the ratio 


rT 


Tqcd 

Aqcd 


^ 0 . 2 . 


(4.4) 


If we dehne low temperature in relation to the critical temperature, let us say, T < 
QATqcd, in terms of the QCD scale, low temperature means 


T < 0.08 Aqcd ~ 0.08 • dvrF. 


(4.5) 
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Returning to our systems in d=2+l with a spontaneously broken rotation symmetry, 
we conclude that the low-temperature region where the effective expansion applies, is 
given by 


T < 0.08 Jss 0.08-21 


(4.6) 


In order to check consistency, consider the d=2+l quantum XY model that is char¬ 
acterized by the Kosterlitz-Thouless phase transition which marks the breakdown of 
the spin-wave picture. On the square lattice, the Kosterlitz-Thouless transition occurs 
at Tkt = 0.343J 97|]. If we define low temperature again in relation to the critical 
temperature, T < QATkt-, in terms of the microscopic scale J, we get T < 0.1 J 
which is perfectly consistent with Eq. 04.61) . 


Regarding the magnitude of the Goldstone boson mass (and hence the external 
held Hs), recall that in QCD the pion masses are roughly Mn ~ 140MeK, leading to 
the ratio 

rM = -/^~0.1. (4.7) 

In fact, the effective theory of QCD even works for kaon masses {Mk ~ SOOMeR) 
and the r^-mass {Mr^ ^ 550MeR) - these particles, much like the pions, represent 
Goldstone bosons. We thus conclude that the effective series derived in the present 
paper are certainly valid in the region 


H„M^ < O.IJ ss 0.1 ■ 2^rF^ 


(4.8) 


which dehnes weak external held Hg, or low mass M.^. Again, the two quantities 
are connected by Eq. 03.3p . Note that the Goldstone boson mass is given in units of 
[Mt^v^] = J, where the spin-wave velocity v has been set to one. 


However, there is an important caveat: although the ratio Hg/J must be small, 
it cannot be arbitrarily small, because of the subtleties arising in two spatial dimen¬ 
sions. In the absence of an external held, spontaneous symmetry breaking does not 
take place in two spatial dimensions at hnite temperature (Mermin-Wagner theorem 
(^b Formally this means that the limit 77^ —)■ 0 in the low-temperature series is 
problematic, as we have discussed on various occasions [3, [sj : if the external held is 
very weak, we leave the domain of validity of the present ehective expansion. Still, 
the excluded region in the plane Hg/J versus T/J is tiny, as we illustrate in Fig. 0] 
for the square lattice quantum XY model and the Heisenberg antiferromagnet. In the 
parameter region above these curves, the ehective analysis presented here perfectly 
applies. In Sec. 14.31 where we consider the staggered magnetization, we explain how 
we have generated these curves. 


In order to illustrate the restrictions imposed by very small Goldstone boson 
masses for general systems (not necessarily spin models where the is determined 


by the staggered held), we follow another strategy. According to Ref. [99|, llQ0|, for 


N > 3 a. nonperturbative mass gap m is generated for the Goldstone bosons at hnite 
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Figure 4: Domain of validity of the effective low-temperature expansion of the square 
lattice quantum XY model (bottom) and the square lattice Heisenberg antiferromag- 
net (top). The effective expansion is safe in the parameter region above these curves. 


temperature for d=2+l dimensional 0(iV) symmetric theories, 


8 \^ T /27rXF2 
1 ) r(l + X) V T 


X 

exp 


2'KXF^^ \ 


1 T 

2 27rF2 



which depends on N through the parameter X, 


X 


1 

X-2 ■ 


(4.10) 


As discussed in detail in Ref. [^, the correlation length of the Goldstone bosons 
^ = 1/Mt^ must be much smaller than the correlation length ^np = 1/'^ defined 
by the nonperturbatively generated mass, in order for the effective expansion to be 
valid. This condition then leads to a relation between the two ratios T/2nF'^ and 
M7r/27rF^, allowing one to identify the parameter domain where the effective series 
are valid. This is illustrated for various values of N in Fig. [5l The region below these 
curves corresponds to the parameter domain where the effective analysis presented 
here cannot be applied. While this excluded region is tiny for N= {3,4,5,8}, it grows 
as N increases. 


We emphasize that the spontaneously broken internal symmetry 0(X) —)• O(X-l) 
not necessarily refers to quantum spin systems. Just for concreteness we continue 
using magnetic terminology, assuming that the underlying system is described by the 
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Figure 5: Domain of validity of the effective low-temperature expansion of d=2+l 
systems with a spontaneously broken symmetry O(A^) —)■ O(A^-l). The effective ex¬ 
pansion is safe in the parameter region above the curves that refer to iV= {3,4,5,8}, 
from bottom to top in the hgure. 


d=2+l antiferromagnetic Heisenberg Hamiltonian for spins in arbitrary dimension N, 
and that the underlying scale can be identihed with the exchange integral J. 


4.2 Manifestation of the Interaction in the Pressure 


If the interaction contribution in the pressure has a positive (negative) sign, we are 
dealing with a repulsive (attractive) interaction. In terms of the dimensionless ratio 


T 

A4’ 

the pressure, up to three-loop order and for general N [N > 2), reads 


(4.11) 


P = 


+ 

+ 


iV- 1 


ho{(j)T^ - 


{N-l){N-3) 


2rTiA 


2 8F2r2 

(AT - l)(iV-3)(5iV- 11) 


hiiaYT 




1287rF4r3 ^ 

{N — 1){N — 3) Pu f 1 I \rp5 I 

-- hi(^) h2[a)T -j^q{a)T + 0{T ). 


(4.12) 
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Figure 6: [Color online] Manifestation of the interaction in the pressure of (2+1)- 
dimensional systems with a spontaneously broken symmetry O(A^) —?• O(A^-l) for 
N= {2,3,4,5,6} from bottom to top in the hgure (vertical cut at mJ/2'kF‘^ = 0.02) in 
the hgure. The curves have been evaluated at T/2ttF‘^ = 0.05. 


As the system is homogeneous, the pressure is obtained from the free energy density, 
Eq. dSU), by 

P = Zo-z, (4.13) 

where zq is the free energy density at zero temperature. 


The low-temperature series for the pressure allows us to extract the information 
on the Goldstone boson interaction we are interested in. The following statements 
apply to (2+l)-dimensional (pseudo-)Lorentz-invariant systems with a spontaneously 
broken internal symmetry O(A^) —)■ O(iV-l). In Fig. [H] we have plotted the sum of the 
two-loop (oc T^) and three-loop (oc T^) contributions to the pressure as a function of 
the ratio Mt^/2ttF‘^. In connection with spin systems, this explicit symmetry breaking 
parameter measures the strength of the external staggered held. Note that we have 
scaled the pressure by the number of Goldstone bosons that appear in a given system, 
i.e., we have plotted the dimensionless quantity 


1 P 
N -IF^' 


(4.14) 


For the temperature T /2ttF‘^ = 0.05 that we have chosen in Fig. [HI the interaction in 
the pressure is repulsive for any of the systems N= {2,3,4,5,6} considered. 


However, the situation becomes more subtle if we go to lower temperatures. For 
A^={2,3}, the interaction remains repulsive, independently of the ratio T/2tiF‘^. But 
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2nF^ 

Figure 7: Manifestation of the interaction in the pressure of (2+l)-dimensional sys¬ 
tems with a spontaneously broken symmetry 0(4) —)• 0(3). The curves have been 
evaluated at the temperatures T/27rF^ = {0.028, 0.025, 0.022, 0.019, 0.016}, from 
top to bottom in the figure. 


the behavior of the systems with > 4 is more complicated as we illustrate in Fig. [7] 
for the case iV=4. By reducing the temperature, the Goldstone boson interaction 
in the pressure becomes attractive in an intermediate region that depends on the 
parameter This region keeps growing as the temperature drops. 

These effects are an immediate consequence of the A^-dependence of the low- 
temperature expansion for the pressure. At the two-loop level, there is just one term 
(oc T^) in Eq. fl4.12p . As it involves a factor of A^-3, one may identify three classes of 
systems regarding the Goldstone boson interaction: the Heisenberg antiferromagnet 
{N=3) where the interaction term vanishes, the quantum XY model {N=2), where 
the interaction is repulsive, and systems with N > 4 where the interaction in the pres¬ 
sure is attractive at two-loop order. Then, at the three-loop level, the Y-dependence 
is more complicated, but overall leads to positive contributions to the pressure. The 
Heisenberg antiferromagnet receives a repulsive contribution, while the repulsive in¬ 
teraction in the quantum XY model gets enhanced. For systems with Y > 4, the 
three-loop corrections are quite large at low temperatures, and involve opposite signs 
compared to the two-loop contribution. This leads to the interesting behavior we 
discussed before, and which we have depicted again in Fig. [8] for the temperature 
T/2ttF‘^ = 0.02 for Y= {2,3,4,5,6}. For quantum spin systems with Y > 4, depend¬ 
ing on the strength of the external staggered held and the temperature, the Goldstone 
boson the interaction in the pressure may be attractive, repulsive, or zero. 
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Figure 8: [Color online] Manifestation of the interaction in the pressure of (2+1)- 
dimensional systems with a spontaneously broken symmetry O(A^) —)■ O(A^-l) for 
N= {2,3,4,5,6} from top to bottom in the hgure (vertical cut at Mt,/2'kF‘^ = 0.02). 
The curves have been evaluated at T/27iF^ = 0.02. 

Note that the curves in Fig.[7]and Fig. |8]do not diverge in the limit Mt,/27iF^ —)■ 0, 
i.e., when the external held is switched off: they tend to a hnite value as in Fig. [6] 
where the temperature is more elevated. This value is positive, except for N=2, where 
the Goldstone boson interaction tends to zero and the system behaves as ideal gas. 
While the limit Mt^/2'kF’^ —)■ 0 is mathematically well-dehned in the series for the 
pressure, one should keep in mind that, due to the Mermin-Wagner theorem, or the 
nonperturbatively generated mass gap, one eventually leaves the domain of validity of 
the effective expansion if decreases. However, for the temperature T/2ttF‘^ = 0.02 
chosen in Fig. [HI the parameter Mt^/2'kF'^ must be really tiny in order for the effective 
series to break down. This can be easily verihed by inspecting Fig. O Therefore, 
although it is conceptually inconsistent to take the limit M^/27rF^ —)■ 0, numerically 
the pressure is not affected. 

We thus emphasize that the effects described above occur in a parameter region 
where our approach perfectly applies. The systematic effective held theory analysis 
predicts that for systems with > 4, depending on the ratios T/2ttF‘^ and Mt^/2ttF'^, 
the Goldstone boson interaction in the pressure may be repulsive, attractive or zero. 
On the other hand, for N= {2,3}, the interaction always is repulsive. 

Finally we mention a subtlety related to the quantities M and TG that both 
depend on the external held Hg, and are related by Eq. fl3.3p . The point is that 
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the interaction already manifests itself at zero temperature: in the presence of an 
external held, the bare Goldstone boson mass M = y/T^^Hl/F receives one-loop and 
two-loop contributions: M —)■ M^. In the hgures shown in the present subsection, 
we have plotted the finite-temperature interaction contribution in the pressure as a 
function of the quantity M^/27rF^ that involves M^, and thus takes into account the 
effect of the interaction at zero temperature. The hgures for the pressure, that all 
refer to hnite temperature, have thus to be interpreted as follows: we start at T=0 
and switch on the external held. Then we switch on the temperature (e.g., in Fig. [S] 
the temperature is T/2tiF‘^ = 0.02), and observe whether the interaction at finite 
temperature is repulsive, attractive, or zero. 

Note again that the systems under consideration not necessarily are quantum spin 
models - the ehects discussed above apply to any (pseudo-)Lorentz-invariant system 
in two spatial dimensions exhibiting a spontaneously broken internal symmetry 0(iV) 

0{N-1). We hnd it remarkable that the ehects concerning the nature of the 
interaction follow from symmetry considerations alone: they neither depend on the 
concrete realization of the system nor on its microscopic details. 


4.3 Manifestation of the Interaction in the Staggered Mag¬ 
netization 


The staggered magnetization at hnite temperature is dehned by 




dz 


(4.15) 


It represents the order parameter, signaling that the internal rotation symmetry is 
spontaneously broken0 


Using the expression for the free energy density fl3.ll) . the low-temperature expan- 

^Strictly speaking, there is no spontaneous symmetry breaking at finite temperature and zero 
external field in two spatial dimensions due to the theorem by Mermin-Wagner. But the staggered 
magnetization still takes nonzero values at finite and low temperatures and weak external field. In 
this sense we refer to the quantity Ss(T, Hs) as order parameter. 
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sion takes the form 


Es{T,Hs) = S,(0,g,)- ^'^ ^hi{a)T 


2 F2 

{N-l){N-3) E ,6 
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hi{af - ^hi{a)h2{a) \T'^ 


(N - l)iN - ?,)ibN - ll)T.sh ( 3 , , ,2 2 , , ,, , .'i „ 

(iV-l)(Ar-3)(3Ar-7) e,6 




+ 


48 F 6 

(iV-l)(Ar-3)2s,6j' 2 

16 

s.t a9(<^)j,3 + o,r4)_ 


|fti((T)’ - -^h^(<Tf h2{a)^T^ 


hi(cr) h2(a) - —hi{a)h2{aY - -;^hi{a) h^ia) [>T'* 


Sir'^F^a da 
where b is given by 

dM"^ 


(4.16) 


167rF3 
b 

ko = 2(/c2 — ki) + — + 


3{N - 3)^/E:^ ^ ^ 


i?4 


F2 647r2F2' 


(4.17) 


The next-to-leading order effective constants ki and ^2 are dehned in Eq. fl2.3p . while 
the emlicit expressions for the coefficients bi and 62 can be found in Appendix B of 

Ref. 0. 

In the two hgures of this subsection, we plot the sum of the two- and three-loop 
contributions in the dimensionless quantity 


= 


. E,(T,hf,)-E,(0,i/,) 

A^ - 1 S. 


(4.18) 


which measures the strength and the sign of the Goldstone boson interaction in 
the temperature-dependent part of the staggered magnetization. Note that the T- 
independent part of the staggered magnetization in the presence of an external held 
is given by 


E,(0,if,) = S, + s/Ws + 0{Hs), 

S, = S,(0,0). (4.19) 

Since we are at zero temperature, the limit —)■ 0 in Eq. fl4.19p raises no problems: 

the Mermin-Wagner theorem does not apply. 

As can be appreciated in Fig. [91 the d=2+l quantum XY model (or any other 
physical realization of N=2) is quite different from systems with N > 3: while the 
sign of the interaction is negative in most of the parameter domain, for small external 
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Figure 9: [Color online] Manifestation of the temperature-dependent interaction con¬ 
tribution in the order parameter (staggered magnetization) of (2-|-l)-dimensional sys¬ 
tems with a spontaneously broken symmetry 0(iV) —)■ O(A^-l) for N= {2,3,4,5,6} 
from bottom to top in the figure (vertical cut at = 0.08). The curves have 

been evaluated at T/2ttF'^ = 0.05. 

helds (small M^) it is positive. For all other systems {N > 3) the temperature- 
dependent interaction contribution in the order parameter tends to large negative 
values if the external held is weak. 

Note that Fig. IH] corresponds to T/27iF‘^ = 0.05. If we go to lower temperatures, 
effects similar to those observed for the pressure take place: while the properties of 
the quantum XY model and the Heisenberg antiferromagnet qualitatively remain the 
same, systems with N > 4 develop positive values for the interaction part in the 
staggered magnetization, as we illustrate in Fig. [10] that refers to T/2tiF‘^ = 0.02. 
In particular, since the curves for N > 4 now exhibit a maximum, the staggered 
susceptibility for these systems vanishes at this point in parameter space (see next 
subsection). Remarkably, in the case of the Heisenberg antiferromagnet, the quantity 
S'Es always takes negative values and never presents an extremum, independently of 
the temperature or the external held. 

The sub leading ehects discussed here are quite small, and the behavior of the 
staggered magnetization is in fact dominated by the one-loop contribution. As one 
would expect, the one-loop contribution in 611 g is negative for any of the systems: the 
corresponding term, proportional to one power of the temperature in Eq. fl4.16p . is 
governed by —hi (which is negative), where hi is the kinematical function that we 
have depicted in Fig. [2] If the temperature is switched on while the external held 
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Figure 10: [Color online] Manifestation of the temperature-dependent interaction con¬ 
tribution in the order parameter (staggered magnetization) of (2-1-1 )-dimensional sys¬ 
tems with a spontaneously broken symmetry 0(iV) —)■ O(A^-l) for N= {2,3,4,5,6} 
from bottom to top in the figure (vertical cut at = 0.03). The curves have 

been evaluated at T/2nF‘^ = 0.02. 

is held hxed (M^ constant), indeed the order parameter decreases. Remember that 
we first switch on the external field Hg at T=0, and then go to finite temperature, 
keeping Hg fixed. 

Interestingly, the effect of the Goldstone boson interaction, as witnessed by Fig. [9] 
and Fig. [TOl does not necessarily follow this behavior observed for the one-loop con¬ 
tribution. For the d=2+l quantum XY model at low temperatures and weak external 
held, the effect is just the opposite: the interaction at hnite temperature leads to an 
increase of the staggered magnetization, meaning that here the interaction tends to 
reinforce the antialignment of the spins and enhances the staggered spin pattern. Sim¬ 
ilarly, for systems with iV > 4, the interaction leads to a slight increase of the order 
parameter at low temperatures for not too weak external helds (see Fig. [TU]) . Overall, 
of course, the one-loop contribution dominates, such that the order parameter still 
decreases when the temperature is switched on. 

It should be pointed out that the low-temperature series for the staggered magne¬ 
tization - in contrast to the series for the pressure - formally diverges if the external 
field is switched off. As we have explained before, taking this limit in any observable 
is problematic, as we leave the domain of validity of the effective expansion. Tak¬ 
ing this limit in the pressure did not numerically affect the series. However, in the 
staggered magnetization and the staggered susceptibility (see below), this limit is di- 
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vergent. But note that the curves displayed in Fig. [9] and Fig. [TO] can be trusted to 
very low values of M^/27rF^ (very weak external held) according to Fig. 0, such that 
the hnite-temperature interaction contribution in the order parameter indeed tends 
to large positive (negative) values for N=2 {N > 3) in a regime where the effective 
expansion still applies. 

At the end of this section, let us explain how we have generated the curves con¬ 
cerning the range of validity of the effective expansion that we have depicted in Fig. jU 
For a hxed value of the external held Hg, the order parameter ^^(T, Hg) drops to zero 
when the temperature is raised sufficiently: at this point the ehective expansion cer¬ 
tainly no longer applies - we are now in a region of elevated temperature where the 
Goldstone boson picture ceases to be valid. To estimate this ’’critical” temperature 
for a given constant value of Hg, we take the hrst two terms in the expansion for the 
staggered magnetization, Eq. fl4.16p . and solve the equation 

TVr 1 AT 1 /-» 

E,(T, //,) = E, + -jrVH.- r ^ 0, (4.20) 

with the kinematical function hi 

hi{a) = —^In fl - ^ —^In fl - . (4.21) 

For a given ratio of Hg/2'KF‘^ (i.e., Hg/J) we thus obtain a ’’critical” temperature 
Tc/J that signals the breakdown of the ehective expansion. While in Fig. IHwe have 
shown the situation for the d=2-|-l quantum XY model and the Heisenberg antiferro- 
magnet (where the specihc values for and the spin-wave velocity Vg have been 

measured), the situation for iV > 4 is similar: the ehective expansion only breaks 
down when the external held is tiny, as depicted in Fig. O 


4.4 Manifestation of the Interaction in the Staggered Sus¬ 
ceptibility 


In connection with spin systems, the staggered susceptibility 


xiT.Hg) 


dT.g{T,Hg) 

dHg 


E, aS,(T,M) 
F2 


(4.22) 


is the response of the staggered magnetization to the applied external held. Rather 
than providing lengthy expressions for xiT, Hg) that can trivially be obtained from 
the low-temperature expansion of the order parameter, Eq. fl4.16p . we just exhibit the 
the general structure of the series: 


x(r, Hg) = x(0, Hg) + xi(r)^ + X 2 (r) + Xz{r)T + 0{T ^). (4.23) 
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Figure 11: [Color online] Manifestation of the temperature-dependent interaction con¬ 
tribution in the staggered susceptibility of (2-1-1 )-dimensional systems with a sponta¬ 
neously broken symmetry 0(iV) —)■ O(iV-l) for N= {2,3,4,5,6} from bottom to top 
in the hgure (vertical cut at M^/27rF^ = 0.06). The curves have been evaluated at 
T/27rF2 = 0.05. 


The coefficients y*, much like the analogous coefficients in the pressure and the stag¬ 
gered magnetization, depend in a nontrivial way on the dimensionless ratio r = 

Regarding quantum spin systems, the staggered held Hg = (0,..., 0, Hg) tends to 
align the spins in an antiparallel pattern and hence increases the staggered magnetiza¬ 
tion. One thus expects the staggered susceptibility y(T, Hg) in general to be positive. 
This is indeed true for the noninteracting (one-loop) contribution in Eq. 04.231) . i.e., 
the leading temperature-dependent term that involves an inverse power of T. How¬ 
ever, the effects related to the interaction are more subtle. 


In the subsequent two hgures, we plot the sum of the two- and three-loop contri¬ 
butions in the dimensionless quantity 


^ TF^ x{T,Hg)-x{0,Hg) 
N-1 S2 


(4.24) 


which measures the strength and the sign of the Goldstone boson interaction in the 
temperature-dependent part of the staggered susceptibility. 


The behavior of the d=2+l quantum XY model is peculiar, as the hnite-temperature 
interaction contribution in the staggered susceptibility - in contrast to all other sys- 
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Figure 12: [Color online] Manifestation of the temperature-dependent interaction con¬ 
tribution in the staggered susceptibility of (2-1-1 )-dimensional systems with a sponta¬ 
neously broken symmetry 0(iV) —)■ O(iV-l) for N= {2,3,4,5,6} from top to bottom 
in the hgure (vertical cut at M^/27rF^ = 0.04). The curves have been evaluated at 
T/27rF2 = 0.02. 

terns {N > 3) - becomes negative for small values of Mt^/27iF‘^, as shown in Fig. [TT] 
and Fig. [T2j The d=2+l Heisenberg antiferromagnet also plays a special role as it 
is the only system with a positive interaction contribution in the staggered suscepti¬ 
bility, independent of the strength of the staggered held (or magnitude of M,r/2vrF^) 
and temperature. For iV > 4 the behavior at T/2ttF‘^ = 0.05 is qualitatively the same 
as for the antiferromagnet (see Fig. m, but for lower temperatures the quantity 6x 
becomes negative for not too weak external helds (see Fig. [12]). But note that for 
small values of Mt^/27iF‘^, the staggered susceptibility (iV > 3) tends to large positive 
values for any temperature - the d=2+l quantum XY model really is the exception. 

Analogous to the staggered magnetization, the manifestation of the Goldstone 
boson interaction in the staggered susceptibility is thus quite subtle. In many cases 
the temperature-dependent interaction contribution in x(T, Hg) is positive, implying 
that the external held reinforces the staggered spin pattern. In some cases, however, 
the T-dependent interaction part in x(T, Hg) is negative, meaning that here the ehect 
of the staggered held is to perturb the antialignment of the spins and to decrease the 
staggered magnetization - or the order parameter in general. This happens for the 
d=2+l quantum XY model at weak external held (and hnite temperature), and for 
systems with Y > 4 at low temperatures and not too weak external helds. It should 
be pointed out again that these subtle ehects regarding the nature of the interaction 
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at finite temperature are rather weak - the noninteracting part in fact dominates, and 
leads to positive staggered susceptibilities, as one would expect. 


5 Conclusions 


Our main results concern the universal nature of the Goldstone boson interaction 
in d=2+l (pseudo-)Lorentz-invariant systems with a spontaneously broken internal 
rotation symmetry 0{N) —)■ O(A^-l). Our effective held theory calculations are valid 
at low temperatures and in presence of a weak external held that explicitly breaks 
internal rotation invariance and thus invalidates the Mermin-Wagner theorem. As 
concrete physical realizations we have considered quantum spin systems. 

The d=2+l quantum XY model and Heisenberg antiferromagnet - or any other 
physical realizations oi N = {2, 3} - are rather special, since the interaction in the 
pressure is always repulsive, independent of the ratio T/2nF‘^ and Mn-/27rF^ that 
measure the temperature and the strength of the external held, respectively. For all 
other systems {N > 4), the interaction in the pressure may be repulsive, attractive, or 
zero, depending on the strength of the external held and the temperature that hence 
control the nature of the Goldstone boson interaction. 

Regarding the staggered magnetization ^^(T, Hg) - the order parameter - the in¬ 
teraction may lead to interesting and subtle ehects. When the temperature is switched 
on while keeping the strength of the external held hxed, in most scenarios the inter¬ 
action tends to decrease the staggered magnetization Ss(T, Hg), as one would expect. 
However, the d=2+l quantum XY model represents an exception, because in very 
weak external helds and at hnite temperature, the interaction tends to reinforce the 
staggered pattern, rather than to destroy the antialignment of the spins. Systems 
with Y > 4 exhibit a similar behavior at very low temperatures in case the external 
held is not too weak: the order parameter is increased in this parameter regime, when 
we go from zero to hnite temperature while keeping the staggered held constant. 

These ehects also manifest themselves in the staggered susceptibility x(T, Hg). In 
weak external helds, the temperature-dependent interaction part in x(T, Hg) tends 
to positive values for N > 3, but in the case of the d=2+l quantum XY model it 
becomes negative. For systems with Y > 4, the T-dependent interaction contribution 
in the staggered susceptibility takes positive values, but for very low temperatures it 
may drop to negative values provided that the external held is not too weak. 

Although the calculation was performed within a (pseudo-)Lorentz-invariant frame¬ 
work, it readily applies to the d=2+l quantum XY model or the d=2+l Heisenberg 
antiferromagnet. The crucial point is that the ehective theory predictions regard¬ 
ing the nature of the Goldstone boson interaction are not ahected by Lorentz-non- 
invariant terms, that would indeed be present in the next-to-leading order Lagrangian 
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and that have not been taken into acconnt. These terms are only relevant for a 
one-loop graph and are contained in the temperature-independent quantity Only 
in this trivial way, and only at next-to-next-to-leading order, do these anisotropies 
emerge in our expressions. In particular, they cannot change the sign of the interac¬ 
tion - and hence alter our conclusions ~ in any of the observables considered. In this 
sense, the effective held theory predictions are universal or model-independent. We 
hnd it indeed quite remarkable that the nature of the Goldstone boson interaction is 
fully determined by the symmetries of the system. 

While we have explicitly referred to the quantum XY model and the Heisenberg 
antiferromagnet - or spin systems characterized by higher values N of spin vector 
components - the results presented here go beyond quantum spin models. Our low- 
temperature series are valid for any d=2+l (pseudo-)Lorentz-invariant system with 
a spontaneously broken internal rotation symmetry 0(iV) —)■ O(iV-l) in presence of 
a weak external held. It is our hope that these subtle ehects may be conhrmed in 
future studies, e.g., by Monte Carlo simulations. 
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